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PROBLEM 35 
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r ; THE ULTIMATE NUMBER @. 


Nine numbers are represented on the cover of this 
issue. Besides the familiar Pi and e, there is the 
Golden Mean; (1 +*/5)/2 and Euler's constant: 


lim (1+ 1/2 + 1/3 + 1/4 + ... + 1/n - In n) 


n—-oc 


plus the DIS numbers P2, X2, F, and G, described in PC8-14. 
The DIS number PR6 is the DIS sum of prime numbers for 
which p+6 is also a prime. 


The first eight of these numbers are known to 30 
significant digits or more; the number PR6 is not presently 
Known at all. The nine components between them include 
representatives of the interesting numbers: irrational 
and transcendental numbers, power functions, the Fibonacci 
sequence, factorials, and primes. 


We seek the product of these nine numbers and make 
the following offer: a year's subscription to POPULAR 
COMPUTING for each new significant digit of the Ultimate 
Number found and verified, up to a maximum of 20. ( The 
limitation is made in case someone is able to calculate the @ 
DIS sum for PR6 analytically. We think that it will 
probably be obtained by actual summation of the primes, and 
hence that the digits of the Ultimate Number will emerge 
one by one.) A digit will be considered verified when 
the subsequent digit is found. 


To clarify the DIS process as it applies to-the 
number PR6, For the l-digit primes, we have only 


205 
-O7 for a sum of .12 


For the 2-digit primes, there are 20 that contribute to 
the sum: 
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0097 for a sum of .0972 


There are 98 3-digit primes that contribute; 584 4-digit 
primes; 3660 5-digit primes; and so on. For all the 
primes of the form p+6 through the prime 9973, the DIS 
sum is around .297, Indicating that the final sum may 
converge to around .33 or so. Convergence is assured, 
so the calculation of UN rests on finding PR6. 


DESK CALCULATOR REVIEVV 
Hewlett-Packard HP-45 


With the introduction of the HP-45, the price on 
the HP-35 was dropped to $300. The 45 sells for $400. 
(Thus the rating for the 35, using the scale given in PC-10, 
rises from 4.911 to 6.576.) The 45 is the upgraded 
version of the 35, and hence this review deals only with 
the added features of the 45 (see the review of the 35 
in PC1-6). 


All the functions of the 35 are retained, plus the 
following: factorial (up to 69!); trigonometric functions 
in degrees, radians, or grads; degrees-minutes-seconds to 
decimal degrees, and the inverse; 9 addressable words of 
storage in addition to the 4-level stack storage; metric 
conversion for volume, weight, and liquid measure; polar 
to rectangular coordinate conversion and the inverse; 
percent and change in percent; summation and means and 
standard deviations. The algorithms used are apparently 
the same as the corresponding functions on the 35 or 80; 
for example, sine 4444230 is .4999947727 on both the 
35 and 45. ‘The standard deviation formula uses N(N-1) 
in the denominator. 


The rating scale comes out 9.100 for the HP-45, 
which indicates that the machine offers a great deal of 
calculating power for its price. A well-written manual 
comes with the machine, although it is reduced to fit in 
the machine's box and is thus difficult to read. 
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ART. COMPUTING 2 


“TESTING 


When a computer program is first submitted to the 
machine for assembly or compilation, it is usually in the 
process of debugging (unless, of course, it contains no 
bugs, but that condition is rare). If the assembler or 
compiler does its work well, all the mechanical bugs will 
be revealed. These include such things as: 


. Undefined symbols. 

. Ambiguous symbols. 

. Mixed mode expressions. 

. Lillegal expressions. 

. Impossible paths (that is, statements which 
cannot be reached in execution). 

. Lllegal values. 

Attempted entry into the middle of loops. 


Any of these errors, and others, should inhibit 
execution and cause the printing of a diagnostic message. 
When all such bugs have been corrected, a run can be made 
that will execute, and the program can be said to be the 
solution to some problem. The question is, is it the 
problem that was supposed to be solved? 


It is at that point that program testing begins. 
We must verify not only that our program solves the given 
problem--under controlled conditions--but that it will 
continue to solve that problem as the input data changes. 
This task is largely art and is not trivial. It is the 
task that demands all the ingenuity and low cunning of 
which the programmer is capable; it is the place where 
the programmer can fully demonstrate his superiority to 
the machine. 


The two stages--debugging and testing--overlap, 
since errors in logic are frequently revealed during the 
debugging stage. Nevertheless, the two stages should be 
thought of as distinct, sinee they require radically 
different techniques. The confusion caused by inter- 
mixing the two (which is true of nearly all textbooks that 
consider either stage) is probably responsible for more 
bad computing than any other single 111 in our field. 

The trouble is compounded by the unfortunate fact that 

the concept of program testing is the single most difficult 
notion to teach in all of computing. All beginners (and 
many “experienced" programmers) fall into traps like the 
following: 


1. I wrote the program myself, and I have great 
confidence in me, and so I know it works. 


2. I @id test my program, and it worked (for one 
case) and E'm sure it will work for every case. 


3. Computers don't make mistakes, and my progrem 
prints all its headings correctly, so I know the results 
are correct. 


Any intelligent person is inclined to laugh at such 
naivete (and admittedly the excuses have been worded as 
baldly as possible), but what other reasons can there be 
for the failure to apply even the crudest of test procedures 
to production programs? Everyone reads of the dramatic 
glitches: the IRS bill for $.01; the rebate check from 
the gas company for $999999999; the book club subscriber 
who receives 3000 copies of the book he rejected--this 
list is endless. Who tested those programs? Everyone 
hears about the big, egregious blunders, like the $18 
million aborted missile, but the minor lapses in testing 
simply cause annoyance, followed, sometimes, by efforts 
at reprogramming. Those who interface daily with program 
packages wonder sometimes whether anyone tests anything. 

Or canes. 


Two trivial examples will illustrate what this is 
all about: 


(A) A subroutine is to be written to move 
the words in block A of storage to block 
B and reverse their order at the same time. 


(B) A program is to be written to number all 
roper fractions that are in lowest terms 
a flowchart for this task was shown in 

issue number 3). 


These two tasks call for different testing techniques. 
In task A, the data can be controlled, and predictable 
and predicted results can be listed. In task B, there 
is no data, and the results from a debugged program must 
be compared to an independent calculation (say, by hand) 
in order to validate the program. 


There is given below a list of guidelines for 
program testing. 


1. You should test the program being tested, and 
not some other program that you prefer to test. The 
commonest student error (other than the error of refusing 
to do any testing at ail) is to change the range of the 
problem so that his work is easier; but then what has he 
tested? If the real program with the extended range then 
blows up in his face, who is to blame? 


For example, suppose a program is written to 
generate prime numbers ending in 7 and add their 
reciprocals. Such a program can only be tested by 

trying it in a limited range, where the results have been 
calculated by some other method, and trusting that it will 
function properly in other ranges. 
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3, The results of a test run should be predictable 
and be predicted. That is, the tester should know what 
he expects to see before it is produced. Ghecking results 
after they are produced is conducive to believing what 
one sees, and computer results always look convincing. 


4, The human eye cannot be trusted to check more 
than a few numbers. Any test procedure that calls for 

printing out block @ (25,000 words) to "see if it is now 
in order" will not work. 


5. No program can ever be guaranteed 100%. ine 
goal is to reach a level of confidence at which the person 
testing can sign off the program and stand behind it (with, 
perhaps, disclaimers about ranges he has not explored 


6. he test procedure should be logically independent 
of the program itself. We seek tests that are independent 
checks, in the same sense that the sum of the interior angles 
of a regular polygon must be (N-2)77 radians. 


7. If any part of the program's logic can be 
deleted and the test procedure still works, then it is 
an inadequate procedure. 


8. There must be positive feedback in any test 
procedure; this means that there must be printed (non-zero) 
results from the test runs. 


To repeat: these are guidelines, not rules. One 
can readily devise situations for which some of them would 
have to be relaxed (in particular, the first guideline). 


Consider a textbook example. A program is written 
to count the l-bits in 1000 words of storage. The program 
has been debugged; that is, it executes and prints a 
result. How shall it be tested? It 1s tempting to load 
the 1000 words with all-alike numbers, or with checkerboard 
patterns, but such procedures are weak (although they had 
better work, to be sure). A more satisfying procedure, 

in the sense of raising one's confidence level, would be 

to load the block of storage with 1000 consecutive numbers. 
Now we have the problem of predicting the number of 1-bits 
in 1000 consecutive integers. Since the controlled data 
is systematic, we can expect to find a formula to use to 
calculate the result. This situation is true more often 
than not. When the data is under our control for test 
purposes, systematic data is readily generated and leads 

to results that can be independently calculated. 
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Students Freguentiy ask "If the predicted resulte 
the machine results don't agree, how can you teli whecner 
the trouble lies in the program helng tested or in the _ 
additional code that was written for the test procedure?” 
You can't. But notice that the converse is much neater: 
if the test procedure is intelligently designed and the 
results do agree, then both sets of code have been validated. 


Another example: consider the calculation of 
factorial 10000 {reported in PCH-10). The computer program 
will produce as output a single number of over 35000 digits. 
Before committing the computer time to production, the 
program had to be thoroughly tested. A preliminary run 
was made to produce (1000!) which was validated by compar- 
ison with previous calculations made on a different machine 
by an independent program (see PC2-10). But that was 
hardly sufficient; we were proposing to go from a tested 
result of 2568 digits to one that pushes out into the 
unknown. So, for this project, independent calculations 
were made to determine the exact number of digits to be 
expected; the number of low order zeros; some of the high 
order digits; and some of the low order non-zero digits. 

All of these results checked, so that the final result 
carried a high level of confidence as to its accuracy. 


The most important point about program testing is 
that it should not be a process that begins after debugging. 
Properly constructed programs should have taken into 

account the notion "How will I test this program?" during 
the analysis and flowcharting stages. A tested program 

is reliable, and “reliability is not an add-on feature." 
Anyone who sets out to learn computing should be taught 

the concept of testing from the beginning---from the time 
when they are first shown how to use a computer to sum 

the numbers in a block of 1000 words. 


Errata 


In the definition of the Square Spiral (Problem 24) 
in PC9-2, the neighbors of square #19 should be squares 
6, 5, and 1372 


The attendee list for the 14th Symposium (PC9-6) 
inadvertently omitted the name of Eric Weiss, of Sun 
O11 Company. 
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In a previous issue (PC7-4) a 
Fortran programming problem was present 
with the challenge that it be coded in 
fewer than eleven statementa. The 
program was to accept input consisting 


20 DATA NUMCI 2» LPREVE( 3728672 TVALUE SG 2108558225519 0Ss1/9 
3@ + INPUT /S4HCLS. s BHIG fe TQUT SARC S% 2s 4H! SHAS 4HMOUNS AHT PUs4HRCHA: 
40 it GHSED -8HS sls QHo/5K-24Hs | BHoAHOMOU so AHNT 
5@ + 4H $21 e4H4// 1 sSHBRs 2a 4HOHBRea4HOKENS 4H DOWS4HN AS s 4H 
62 a 48LOWS 2 44<-/ 7) -4G4HOXe Ls 4s iSe4HH DO,AHLLAR;: 4H BIL edHL(S)s 
a’) + 4H/LGXsGH- Leo 2sSH1EK 2s AKHHALF 24H DOL» SHLAR/ Ss 4HI3Ks s4HI 24 Be 
88 + 4HH QU.» 4HARTEs 4HR/13 5 4HX212s4Hs BH 
98 + 4HE227s4HH NI »4HCKELs 4H/1 3X2 4Hs 122 24Hi 1H 
166 + 4HS9 497 
118 58 READ(G@eINPUT) ITENs IPUR 
126 DO 156 1=2,7 
130 EPREVCT=SNUMC 1 <1) *IVALUECI-1)+IPREVCI-1) 
140 158 NUMCI=CITEN-IPUR-IPREVC IJ /SIVALUECI) 
158 WRITECL» IOUT) IPURsITENs CNUMCI)» T=257) 
166 GO TO 58 
176 END 
READY 
RUN 
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of the amount purchased by a hypothetical 
customer and the amount tendered by him 
in payment for the purchase. The 

output was to be the change due hin, 
broken down into the minimum number of 
bills and coins possible. 


One reader was able to code the 
roblem in only nine Fortran statements 
busine proper ANSI Fortran) and is to 
be congratulated for his ingenuity; 
he is Clint Erickson, Jr., 23441 Aetna, 
Woodland Hillis, CA 91364, and his 
solution is given on this page. (Keep 
in mind that the problem statement 
ignored run-time efficiency of the code.) 


‘3. oe 
: ay : ; ; woo z 
. ere Two problews this month, one # 
| return to the COBOL. pvoeram in last month's 26 Cid SUEC2E 
issue. and the other a new Fervvan problem. 28 Le4uCars 


4G GRITEC1 2193 Job 

Ss. Tne COBOL problem presented in PCG-13 oa 18 FORMAT C/14s 2% sA4) 

ed three questions concerning the program 62 STOP 
presented. if the PICTURE clauses for 78 END 
A, B, C, and N were changed to read PICTURE 99, 82 SUBROUTINE SUB C1+sK) 
what would your answers to the three 98 T=1+318656 
questions be? 140 K=1/1068090 

112 RETURN 
The Fortran code (a main program and 12g END 


subroutine) presented here was run on a CDC 3170 READY 
{the code is machine dependent, but this does not RUN 
deter from the problem) and produced the line of 
output shown. Why? 
35 DOGS 
Answers to these problems should be sent 
to Speaking of Languages, POPULAR COMPUTING, STOP @ 
Box 272, Calabasas, CA 91302. 
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The new 208DM, Desk-Top Calculator. 
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*Combination 908DM and Programmer Kit... 
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A COMPUTER PERSPECTIVE 
by the office of Charles & Ray Eames 
Harvard University Press, 1973, 174 pages, $15 


"A Computer Perspective is an illustrated essay on 
the origins and first lines of development of the computer. 

"Me book is based on an exhibition conceived and 
assembled for International Business Machines Corporation. 
Like the exhibition, it is not a history in the narrow sense 
of a chronology of concepts and devices. Yet these pages 
actually display more true history (in relation to the 
computer) than many more conventional presentations of the 
development of science and technology. 

"Te exhibition...was designed to illuminate the 
creative forces and social impact of history being made in 
our own times." 


This is the compubtst ls picture book. it begins with | 
a picture of Babbage's punched card, and the schematic he 
drew of his "Great Calculating Engine.' There are photos 
of every forerunner of the computer, back for nearly 200 
years. 


Here is Herman Hollerith, his original keypunch and 
tabulator, and the patent on the latter. The first Russian 
punched card; early adding machines; a 4-bank set of 
Triumphator calculators; the Powers line of punched card 
machines; the first IBM 100 Percent Club; all the people 
involved with early attempts at automatic computation-- 
this is the scrapbook par excellence. Some of the material 
goes far afield (Landon and Knox campaign buttons; Tik-Tok 
from The Road to Oz) but the authors must have been faced 
with awesome problems of what to include and what to reject. 
Von Neumann's first flowchart is here, as is a flowchart 
of Grace Hopper's of 1949. 


Although the exhibition was commissioned by IBM, the 
photos and exhibits do justice to the contributions of 
Powers, Stibitz, Zuse, Eckert, and Mauchly. The Selective 
Sequence Electronic Calculator gets two full pages, and the 
overall view of it shows the central pillars (which were 


carefully air-brushed out of later photographs). 


The history proceeds up to 1950, and includes ENIAC, 
SEAC, UNIVAC I, the IAS machine, EDSAC, and WHIRLWIND I. 


fhe omissions are interesting. There are no pictures 
of the 6028. 604, or the CPC, each of which was a landmark. 
Software is neglected, except for a reference to COBOL. 
Patents are shown, but not the patent on the rectangular hole. 


After browsing through the book, the reader should then 
browse through its index, to discover more hidden treasures, 
Listings like "Elliot Frog and Switch Company" in the index 
are bound to trigger more browsing. The book is an almost 
endless delight. 


The Eames group has gone to extraordinary pains to 
obtain the original documents and facts about our history. 
This is surely a book that can be recommended as belonging 
on the bookshelf of anyone interested in computing. 
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BIT COUNTING 


PROBLEM 36 


Given a word-oriented binary computer with 32768 
words of core storage. A program is to be written that 
will reside at the top end of core, from word S through 
word 32767. This program will consist of two loops. 
The first loop will load every word of core from 00000 
through S-2 with its own address; that is, word 00001 is 
to contain the number 1; word 00002 is to contain the 
number 2; and so on. The second loop counts the number 
of l-bits in all 32768 words of core storage and puts 
that count in word S-l. 


Problem: To predict the final contents of word S-1. 


The conditions of this problem will vary from 
machine to machine, and with the length of the program. 
The following formula: 


4 OES A 441 4 
at + MAX u(moa it?) - 2 +1, 0| 
ye | pitt 


i=0 


gives the number of l-bits in the integers from 1 to M. 
J is the "width" of Min bits. For example, 


32700), = 111111110111100),, and J = 15. 


The brackets denote "greatest eaters in." For the 
number 32700, the formula yields 244952. 
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120 DEF FNB(M} 


130 FNB=0 
140 G 


160 T= 

170 FOR I=0 ™ J 
180 lect a 
190 T2=FNC(M, 2 


210 T2=0 

220 T=T+T1+Te 
230 NEXT I 
240 FNB=T 
250 ENDFNB 
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TRYING TO FIND GOOD COMPUTER 
LITERATURE. , . 


. AND THE TIME TO READ IT? 


Hire full-time research for just $4.25 a month! Here's what you 


get: 


1. 
2. 
3. 


4. 
5. 


a staff of computer pros continuously monitoring the 
computer literature 


a technical library source of 59 computer publications 
and 123 trade/management publications 

news of conferences, meetings, seminars 

reviews of new books 


original reports about problems faced and solved, but not 
yet reported in the literature 


. «presented in report form cach month. Write for information 
about DATA PROCESSING DIGEST. Or send $4.25 for our 
current issue and apply to your continuing subscription (12 
issucs, $51). 


Data Processing Digest,Inc. 3€ 


6820 LA TIJERA BOULEVARD, LOS ANGELES, CALIFORNIA 00045 / PHONE (213) 776-4334 


The following program will evaluate the bitecount formula 
for the integers from X to ¥. 


100 DEF FNA(X,Y)=FNB(Y)-FNB(X-1) 
120 IF M>O THEN 150 


OTH 250 
150 EEO MEO S) )+1 


atic 
(I+1))-2 
200 IF T2>0 THEN 220 


260 DEF FNC(X, ¥)=X-Y*INT(X/Y) 
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Given a calculating device that can operate on 
8-digit operands to produce an S-digit result. Addition 
and subtraction of larger numbers can be done by parts 
fairly obviously: 


3.1415926 | 53589793 
+2 .7182818 |} 28459045 


5.8598744 | 82048838 


Multiplication of 8-digit numbers to produce a complete 
16-digit product is also quite straightforward: 


Azi41 59268 
62718 28185 


1669 9468 
0885 1338 
1610 6868 


0853 _7238 
853 9733 9875 9468 


If the four parts of the original 8-digit numbers are 
iabelled as shown, then the partial products are, in order, 
BD, AD, CB, and AC, and each of these is 8 digits long, 
leading to the full 16 digit product. 


Extending a quotient can be done as follows. Suppose 
we wish to obtain t digits of the quotient of the numbers 


31415926 
2.7182818 


An 8-digit machine will readily produce the first 8 digits 

of the quotient: 1.1557273. This quotient can be multiplied 
back by the divisor, using the previously explained methods 
of double precision multiplication. The product in this 
case is 3.14159248535314. If this is subtracted from the 
dividend, the remainder is 11464686, and that remainder, 
divided again by 2.7182818, gives the next 8 digits of 

the quotient, 42176222. 


All of which leads to a method of extending a divisor. 
Label the parts of the division as follows: 


322415926 53589793 S A+B 
2.7182818 28459045 % Gees) 
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The calculation of A/C must be carried to 16 places, as 
outlined above. The other two terms need only be carried 
to 8 digits. The final calculation, then, is: 


1.1557273 42176222 
+ 19714583 
- 12099884 


1.1557273 . 49790921 


which is correct to 16S. 


In the early days of computers, when the operation 
code DIVIDE was a costly frill, special methods were 
devised to yield a quotient by iterative schemes. For 
example, the Newton-Raphson scheme can be applied to 
the equation 1/x - N= 0, leading to the recursion: 


Be = x (2 - Nx,,) 


n+1 


where N is a number whose reciprocal is desired. 
expressed in the form 


pi Bee x, (1 - Nx,,) 


then the only double precision calculation required is 
in the term Nx, 


For example, to obtain the reciprocal of N= 7, 
using a starting value of Xo = «1: 


Ky =. + (2b - 7) 5.13 
Ep = .13 + .13(1 - .91) = .1417 
x3 = «1417 + .1417(1 - .9919) = .14284777 


As a method for extending a divisor on a pocket calculator, 
this would be a slow and very tedious process. 
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